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ABSTRACT 



We report on the measurement of the physical properties (rest frame K-band lu- 
minosity and total stellar mass) of the hosts of 89 broad line (type--l) Active Galactic 
Nuclei (AGN) detected in the zCOSMOS survey in the redshift range 1 < z < 2.2. The 
unprecedented multi-wavelength coverage of the survey field allows us to disentangle 
the emission of the host galaxy from that of the nuclear black hole in their Spectral 
Energy Distributions (SED). We derive an estimate of black hole masses through the 
analysis of the broad Mgll emission lines observed in the medium-resolution spectra 
taken with VIMOS/VLT as part of the zCOSMOS project. We found that, as compared 
to the local value, the average black hole to host galaxy mass ratio appears to evolve 
positively with redshift, with a best fit evolution of the form (1 -|- z)^'^^^'^'^'^-"-^ , where 
the large asymmetric systematic errors stem from the uncertainties in the choice of IMF, 
in the calibration of the virial relation used to estimate BH masses and in the mean 
QSO SED adopted. On the other hand, if we consider the observed rest frame K-band 
luminosity, objects tend to be brighter, for a given black hole mass, than those on the 
local Mbu-Mk relation. This fact, together with more indirect evidence from the SED 
fitting itself, suggests that the AGN hosts are likely actively star forming galaxies. A 
thorough analysis of observational biases induced by intrinsic scatter in the scaling re- 
lations reinforces the conclusion that an evolution of the Mbh — relation must ensue 
for actively growing black holes at early times: either its overall normalization, or its 
intrinsic scatter (or both) appear to increase with redshift. This can be interpreted as 
signature of either a more rapid growth of supermassive black holes at high redshift, a 
change of structural properties of AGN hosts at earlier times, or a significant mismatch 
between the typical growth times of nuclear black holes and host galaxies. In any case, 
our results provide important clues on the nature of the early co-evolution of black 
holes and galaxies and challenging tests for models of AGN feedback and self-regulated 
growth of structures. 

Subject headings: galaxies: active - galaxies: evolution - quasars: emission lines - cos- 
mology: observations 
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1. Introduction 

Tight scaling relations between the central black holes mass and various properties of their 



terize the structure of nearby inactive galaxies 


(Maeorrian et al. 
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2009). A result of the search for local OSO rehcs 



via the study of their dynamical influence on the surrounding stars and gas made possible by the 
launch of the Hubble Space Telescope nearly twenty years ago, these correlations have revolution- 
ized the way we conceive the physical link between galaxy and AGN evolution. Coupled with 
the fact that supermassive black holes (SMBH) growth is now known to be due ma inly to radia- 
tively efficient accret i on over cosmological times, taking place durin g "active" phases (jSoltan Ill982l : 
Marconi et al. II2OO4I : Ishankar et al. II2OO4I : iMerloni fc Heinz Iboosi ). this led to the suggestion that 
most, if not all, galaxies went through a phase of nuclear activity in the past, during which a strong 
physical coupling (generally termed 'feedback') must have established a long-lasting link between 
host's and black hole's properties. 

This shift of paradigm has sparked the activity of theoretical modelers. Following pioneering 
analytic works (Ciotti & Ostriker 1997,2001; Silk & Rees 1998; Fabian 1999; Cavaliere & Vittorini 
2002; Whyithe & Loeb 2003), widely different approaches have been taken to study the role of AGN 
in galaxy evolution. Semi-analytic models (SAM) have been the most numerous (Monaco et al. 
2000; Kauffmann & Haehnelt 2000; Volonteri et al. 2003; Granato et al. 2004; Menci et al. 2006; 
Croton et al. 2006; Bower et al. 2006; Malbon et al. 2007; Marulh et al. 2008; Somerville et al. 
2008), and, among other things, have helped establishing the importance of late-time feedback from 
radio-active AGN for the high-mass end of the galaxy mass function. However, firmer conclusions 
on the physical nature of such a feedback mode have been hampered by the large freedom SAM 
have in choosing baryonic physics recipes to implement in their schemes. On the other hand, fully 
hydrodynamic simulations of the cosmological evolution of SMBH have been also performed (Di 
Matteo et al. 2003, 2008; Sijacki et al. 2007; Goldberg & Di Matteo 2008), but their computational 
costs have so far allowed only a limited exploration of sub-grid prescriptions (AGN physical models) 
in relatively small cosmological volumes. A third, hybrid, approach has also been followed, in which 
the results of high-resolution hydrodynamic simulations of galaxy-galaxy mergers with black holes 
(Di Matteo et al. 2005; Springel et al. 2005) have been used to construct a general framework 
for merger-induced AGN feedback, capable of passing numerous observational tests (Hopkins et al. 
2006). 

Almost all the approaches outlined above use the local scaling relations as a constraint to 
the model parameters. As such, these relations have proved themselves unable to unambiguously 
determine the physical nature of the SMBH-galaxy coupling. One obvious way out of this impasse is 



to study their redshift evolution, that different models predict to be different (iGranato et al. 



2004; 



Robertson et al. 


2006: 


Croton 


2006: 


Fontanot et al. 


2006: 


Malbon et al. 


2007: 


Marulli et al. 


2008; 


Hopkins et al. 


2009). 
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From the observational point of view, the current situation is far from being clear. In recent 
years, a number of groups have employed different techniques to try and detect signs of evolution 
in any of the locally observed scaling relations. Most efforts have been devoted to the study of the 
Mbk — (y* relation. Shields et al. (2003) and Salviander et al. (2007) have used narrow nebular 
emission lines ( [OIII] , [Oil] ) excited by the AGN emission in the nuclear region of galaxies as proxies 
for the central velocity dispersion, and compared these to the black hole mass estimated from the 
broad line width of QSOs from z ~ to z 3 (see In both cases, a large scatter has been 
found in the rela tion between Mrh a nd u*, and the results are either in favor of (Salviander et al. 
2007) or against (jShields et al. 1120031 ) a positive evolution of the black hole mass to host dispersion 
ratio. However, as pointed out by Botte et al. (2004) and Greene h Ho (2005), there are a number 
of problems with the underlying assumption that the narrow emission lines are good probes of the 
central gravitational potential, and the systematic uncertainties this method is endowed with are 
large. Komossa and Xu (2007) have also shown that the the AGN for which the [OIII] emission 
line width is far broader that the host galaxy's stellar velocity dispersion o"* tend to show clear sign 
of blue-shift in the narrow emission line (of the order of 100 km/s or larger). Thus, using [OIII] as 
a proxy for would at least require good enough spectral resolution to measure such blue-shifts. 

An alternative path has been followed by Woo et al. (2006), Treu et al. (2007) and Woo 
et al. (2008), who have studied carefully selected samples of moderately bright AGN in narrow 
redshift ranges {z ~ 0.36 and 0.57), where the host's stellar velocity dispersion can be measured 
directly from the absorption lines in high signal-to-noise spectra. They also found evidence of 
(strong) positive evolution of the A/bh to cr* ratio compared to the local value. This method, 
although promising and reliable, is quite inefficient and telescope-time consuming: secure detection 
of spectral absorption features in massive ell ipticals at 1 < S. 2 require hundreds of hours of 



integration time on a 8-meter class telescope (ICimatti et al. 1120081 ) 



Other groups have chosen to try and derive inform ation on the host mass of broad line 



AGN using multi-co lour image decomposition techniques (IJahiike. Kuhlbrodt fc Wisotzki 



Sanchez et al. I l2004l ) or spatially deconvolving optical spectra (ILetawe et al. 1120071 ). Due to the 



2004 



severe surface brightness dimming effects, employing these techniques for redshift QSOs becomes 
increasingly challenging (but see Schramm, Wisotzki and Jahnke 2008; J ahnke et al. 2009; Decarli 



et al . 2009, in prep.), unless gravitationally lensed QSOs are selected (jPeng et al. I l2006bl : IRoss 



20091 ). In all cases, very deep, high resolution optical images [HST) are necessary to reliably dis- 
entangle the nuclear from the host galaxy emission. Statistically, the most significant results have 
been published by Peng et al. (2006b), who, based on a large sample of 51 AGN (both lensed and 
nonlensed) in the range 1 < z < 4.5, observed that the ratio Mbh/-^* increases with look-back 
time, up to a factor ~ 4^^ at the highest redshift probed. 

Also cold/molecular gas motions on large galactic scales have been used to infer the properties 
(total mass in particular) of the hosts of AGN at different redshifts. In the local Universe, HI 21 
cm lines have been used by Ho et al. (2008) to derive total stellar masses of the bulges around 
Seyfert-like AGN. Using instead CO lines to estimate the velocity dispersion in high-redshift QSO 
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hosts, Walter et al. (2004) and Shields et al. (2006) find tentative evidence that the hosts are very 
under-massive compared to their central BHs (i.e., a positive redshift evolution of the Mbh/-^^* 
ratio), a result confirmed by the study of Ho (2007). 

Finally, a completely different approach has been that of trying to follow and compare the 
evolution of global descriptors of the galaxy and SMBH populations, such as mass densities (Mer- 
loni, Rudnick and Di Matteo 2004; Hopkins et al. 2006b; Shankar et al. 2008). Using the simple 
ansatz that total black hole mass density can only increase with time, and requiring that the limits 
imposed by local demographics are not violated, these works showed in general very moderate, if 
any, signs of cosmological evolution of the average black hole to host mass ratio. 

Here we present a new method to tackle the issue of studying black hole-galaxy scaling relations 
at high redshift. Starting from a sample of un-obscured AGN, for which the broad line kinematics 
can be used to infer the central SMBH mass, we take advant age of the unprecedentedly deep multi- 



wavelength coverage of the COSMOS ([Scoville et al. 1120071 ) field and develop a novel SED fitting 



technique that allows us to deco mpose the entire spectr al energy distribution into a nuclear AGN 



and a host galaxy components (|Bongiorno et al. 1120071 ). We show here how, for the majority of 



the objects in our sample, rest frame K-band luminosity and total stellar mass of the host can be 
robustly determined, opening the way to a detailed study of the scaling relations in type-1 AGN 
at 1 < z < 2.2. 

The structure of the paper is the following: we will begin (§ [2]) by introducing our sample, 
before proceeding to a discussion of our SED fitting method in section [3l focusing our attention on 
the measures of the rest-frame K-band luminosity of the AGN hosts (§ 13. ip . on their total stellar 
mass (§ 13. 2p . In section d] we will describe our estimates of the black hole masses obtained by 
studying the properties of the broad MgH emission line, while in § 14.31 we will briefly outline the 
characteristics of our AGN sample in terms of bolometric luminosity, BH mass and accretion rates. 
Section [5] contains the main novel results of our study, namely the analysis of the scaling relation 
for the objects in our sample, as well as a characterization of their observed redshift evolution. 
An important part of our analysis, however, is the assessment of possible observational biases 
responsible for the trend observed, that we carry out in section [H This allows us to reach robust 
conclusions, that we discuss at the end of the paper, in section [71 

Throughout this paper, we use the standard cosmology {Qm = 0.3, = 0.7, with Hq =70 
km s~^ Mpc~^). 



2. The zCOSMOS type 1 AGN sample 



Our AGN sample con sists of the sub-sample of objects in the zGOSMOS bright spectroscopic 
catalog (jLillv et al. 1120071 ) for which one or more broad emission lines have been identified in the 
spectrum. As such, it will be in the following identified as either a Broad Line AGN (BLAGN) or 
a type-1 AGN sample, without introducing any distinction between the two terms. 
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The zCOSMOS bright sample consists, at the times of writing, of 10,644 (medium resolution, 
MR) spectra observed with the VIMOS multi-object spectrograph on ESO-VLT in the COSMOS 
field, selected only on th e basis of their I a r magn itude (Lab <22.5), based on the HST/ ACS imaging 
of the COSMOS field (IKoekemoer et al. I l2007l ) . Fo r a detailed description of the spectroscopic 



survey we refer the reader to iLillv et al. I (|2007l . |2009|) 



Within the zCOSMOS database, we have selected objects with broad emission lines full width 
half maximum (FWHM) larger than 2000 kms~^, a secure threshold for truly broadened lines, as 
compared to our spectral resolution (R ~ 580 for the MR grism, corresponding to ~ 4.8^ and 
~ 520 km/s at the wavelength of Mgll emission). The final sample of type-1 AGN spectra selected 
from the zCOSMOS survey consists of 164 objects which correspond to about 1.8% of the objects 
in the total 10k zCOSMOS database with measured redshift. 

We measure bla ck hole masses by applying the "virial " or "empirically calibrated photo- 



ionization" method (jWandel et al. I Il999l : iKaspi et al. I l200d : iPeterson et al 



2004; Bentz et al. 



20061 ) to the MgIIA2798A broad emission line. We have selected BLAGN in the redshift interval 
z ~ [1.06, 2.19], within which MgIIA27984 can be measured reliably. This takes into account the 
expected broad line width and the problems occurring at the edge of the spectrograph (e.g. fringing 
in the red part of the spectrum) we estimate that the spectral range available to measure Mgll 
line widths is ~ 5650A-9150A (VIMOS spectral wavelengths range from 5500A to 9500A). 

Within this range, the zCOSMOS bright sample contains 104 AGN. After a quick inspection, 
15 of them have been excluded from the analysis because of the low quality of the available spectra, 
leaving us with 89 objects. Ten of those have radio counterparts (at 1.4 GHz) in the VLA/COSMOS 
catalogues (Schinnerer et al. 2007; Bondi et al. 2008), and are listed in Table [2j 



2.1. Photometry 

The zCOSMOS BLAGN sample has been cross correlated with the optical multi-band catalog 
of Capak et al. (2007), the CFHT/K band catalog of Mc Cracken et al. (2009), the IRAC catalog 
by Sanders et al. (2007) and the MIPS catalog by Le Floch et al. (2009). Briefly, the optical 
catalog contains about 3 million objects detected in at least one of the Subaru bands (b,v,g,r,i,z) 
down to a AB magnitude limit of ~ 27 (see Capak et al. 2007, Taniguchi et al. 2007 for more 
details). From this catalog is possible to extract a subsample of ~ 1.3 million sources which have 
signal to noise > 5 in the i- or z-bands. The K-band catalog contains about 5x 10^ galaxies detected 
at a S/N>5 down to K(AB)=23.5 (Mc Cracken et al. 2009). The IRAC catalog contains about 
4x10^ objects detected in the 3.6 micron (IRAC channel 1) band and it is 90% complete at > 1//Jy 
(AB=23.9). For each source in the catalog, the photometry from all the other IRAC channels is 
also reported. The MIPS catalog, obtained in Cycles 2, 3 and 4, has very accurate photometry 
(Sanders et al. 2007). As described in details in Salvato et al. (2009), the fluxes in the optical 
and NIR bands were measured in fixed apertures of 3" diameter, on point-spread function (PSF) 
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matched images (FWHM of l."5). Monte-Carlo simulations (Capak et al. 2007) have been used to 
correct for the flux potentially missed within the apertures. The 3." 8 aperture fluxes given in the 
the COSMOS-IRAC catalog sources were also converted to total fluxes by using conversion factors 
taken from Surace et al. (2005). 

To make sure that confusion is not an issue in the IRAC and MIPS bands, we have visually 
inspected all IRAC and MIPS matches. In the IRAC bands, 82/89 objects have secure counterparts, 
and for them no obvious case of blending has been found when looking into 3" diameter circles 
centered on the source of the optical photometry: only one IRAC source is found in all objects, and 
contaminating flux is usually lower than 1%. For MIPS sources, the situation is slightly different. 
There we found 4 cases where the MIPS source could indeed be a blend of two IRAC sources lying 
within the MIPS error circle. We have thus decided to remove the MIPS photometric point for 
these objects. In any case, we have verified that our results are not changed, neither quantitatively 
nor qualitatively if we recalculate the stellar masses increasing the error uncertainty on all the 
24/im points by it 50%. 

In summary, apart from the 7 cases with ambiguous identification of the IRAC counterparts 
and the 4 confused MIPS sources, for all remaining 78 objects in our final sample, we have used 14 
different bands that encompass optical to MIR wavelengths: 6 SUBARU bands (B, V, g, r, i, z); U, 
J and K band from CFHT + 4 Spitzer/IRAC bands + 24/im form Spitzer/MIPS. This allows us 
to sample a wide wavelength interval, ranging from ~3800 A(Ucfht) to 24^m. All errors quoted 
are Poissonian. 



3. Disentangling the AGN and host galaxy emission with SED fitting 

One of the crucial goals of our study is to use the unprecedented multi-wavelength coverage 
of the COSMOS field to robustly derive host-galaxy properties through detailed model fitting of 
the total SED of the broad line AGN in our sample. Salvato et al. (2009) have demonstrated how 
the COSMOS data allow the determination of reliable photometric redshifts by using composite 
AGN+galaxy templates to fit the multi-band photometry of all XMM-COSMOS sources (including 
both obscured and un-obscured AGN). Here we apply a similar technique to our sample of BLAGN 
(with known spectroscopic redshift) to try to unveil the physical properties of the galaxy component. 
We fit the observed SED with a relatively large grid of models made from a combination of AGN 
and host galaxy templates. For the AGN component we adopt the Richards et al. (2006) mean 
QSO SED (but see ^ 16.31 for a discussion of possible alternative choices), as derived from the study 
of 259 IR selected quasars with both Sloan Digital Sky Survey and Spitzer photome try. We allow 



for e xtinction of the nuclear AGN light applying a SMC-like dust-reddening law (|Prevot et al. 
19841 ) of the form: Ax/E{B — V) = 1.39 A~^-^ for E{B — V) (reddening factor) values in the range 



< E[B — V)< 0.3. For the host galaxy component, we adopt two different sets of templates: 



i) First of all, we use the library of (observed) galaxy templates produced by the SWIRE survey 
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(see Polletta et al. 2007, hereafter P07). From the entire hbrary available of 25 templates, 
we excluded the AGN and composites (starburst+AGN), thus retaining only 14 templates 
(3 Ellipticals, 7 Spirals and 4 Starburst). Such a fitting algorithm has four free parameters: 
two normalizations for the AGN and galaxy templates, respectively, and two corresponding 
reddening factors. 

ii) We also created our own library of synthetic spectra using the well known models of stellar 
population synthesis of Bruzual and Chariot (2003, hereafter BC03). Similarly to previous 
studies of the galaxy population in COSMOS (Ilbert et al. 2009; Bolzonella et al. 2009), we 
built 10 exponentially declining star formation histories (SFH) SFR oc e~*''sc/T" ^j^^j^ e-folding 
times, r, ranging from 0.1 to 30 Gyr, plus a model with constant star formation. For each of 
these SFH, we calculate the synthetic spectrum at different ages, tagej ranging from 50 Myr to 5 
Gyr, subject only to the constraint that the age should be smaller than the age of the Universe 
at the redshift of the source. Finally, we allow for dust extinction, modeled by means of the 
Calzetti's law (Calzetti et al. 2000), with values in the range < E{B — V)< 0.5. Following 
Fontana et al. (2006) and Pozzetti et al. (2007), we impose the prior E{B - V) < 0.15 if 
tagc/T > 4 (a significant extinction is only allowed for galaxies with a high SFR). We adopt 
a Salpeter (1955) Initial Mass Function (IMF) to compute stellar masses. Different choices 
of IMF lead to systematic shifts in t he estimated stel lar masses for any given SED (with the 



maximum shift for a Chabrier IMF ( Chabrier 2002) given by -M^, chabrier 



see e.g. Pozzetti et al. 2007; Ilbert et al. 2009). We will discuss the effect of the IMF on our 
results in section [6l Such a fitting algorithm has six free parameters: two normalizations for 
the AGN and galaxy templates, respectively, two corresponding reddening factors, and the 
age and e-folding time of the exponentially declining SFH. 



We examined the global spectral energy distribution of each object and we fit the observed 
fluxes using a combination of AGN and galaxy emission using the templates extracted from our 
libraries. Nine examples of the SED fitting are shown in Figure [1] (a gallery of all 89 SEDs can be 
found at www.inpe.mpg.de/~ain/plot_sed_all_rev.pdf). 

For each parameter of interest {Mk, Af*), we compute the one-dimensional distribution 
obtained marginalizing over all the other parameters (with flat priors). The normalized probability 
distributions of Mk {P oc exp(— x^d/^)) ^ objects are shown as insets in Fig. [H As a 
general rule, because we are fitting the data with the sum of two model components, the probability 
distributions are asymmetric, with sometimes large tails towards small values of these parameters, 
corresponding to the cases in which the fitting procedure does not require with high significance 
the galaxy component besides the AGN one. We determine the best fit value of the parameter 
of interest as the value that minimizes the x^- One sigma errors on the best flt parameters are 
computed rescaling the observational uncertainties until the minimum reduced Xrcd ~ ^ (f*^^ ^ 
number of degrees of freedom equal to 8 for 79 objects, 7 for 4 without MIPS and 5 for the 
7 without clear IRAC counterparts), and then flnding the range in the parameter values within 
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which Ax^ < 1. We then assign just an upper hmit to the host galaxy rest-frame K-band magnitude 
in all those cases where, within the above-mentioned uncertainty (i.e. where Ax^ < 1), the rest 
frame K-band luminosity can have values smaller than a fixed fraction /gai of the corresponding 
AGN luminosity in the same rest-frame band. By inspection of the best fit SED, after a number of 
trials, we fix /gai = 0.05; simply put, if we find a significant probability that the galaxy component 
is smaller than 5% of the AGN one, we decide that only an upper limit to the galaxy luminosity 
can be assigned. These 5% limit values are marked as vertical red lines in the insets of Fig. [TJ 
Increasing such a threshold to 0.1, although almost doubling the number of objects with only 
upper limits on the host galaxy Mx, hardly produce any significant change in the general results 
and in the global trends discussed in section [5l In those cases, the value of the upper limit is taken 
as the (non-zero) value of the parameter associated to the highest possible value of Mk within the 
uncertainty. Finally, whenever we decide that only an upper limit can be meaningfully associated 
to the K-band luminosity of the host galaxy, we assign an upper limit to the object's total stellar 
mass adopting the median mass-to-light ratio of all other objects in the sample (this corresponds 
to adopting the following relation between the logarithm of the total stellar mass and Mk for the 
upper limits in the sample: LogM,,, = — 0.55 — 0.4(Mx —3.28), see § l3.2p . In total, for 10/89 objects 
we can provide only upper limits for the host galaxy SED component (and thus for Mk and M^). 



3.1. Rest frame K band luminosities 

We are mainly interested here in determining the total mass of the host. It is generally believed 
that local scaling relations apply only when the bulge/spheroid component of the host galaxy is 
considered (Kormendy and Gebhardt 2001; but see the recent works of Kim et al. 2008, Bennert 
et al. 2009, for a different view); however, reliable bulge-disk decomposition for our AGN hosts are 
problematic and will not be considered here; we will briefly discuss the implications of this issue 
later, in section IHTTl In most of the objects of the sample the nuclear AGN emission dominates the 
emission in all optical bands, and the constraints on the host galaxy emission are derived mostly 



in the wavelength ran ge where the AGN SED has a minimum, around 1.2 ^Lim, (jElvis et al. Ill994l : 



Richards et al. 1120061 ). This is close to the rest frame K-band, which is itself a good (i. e. such that 
the mass-to-light ratio in the K band has a l-o" scatter of about 0.1 dex) indicator of total mass (see 
e.g. Madau, Pozzetti & Dickinson 1998, Bell et al. 2003). We thus proceed in two steps: we first 
try and constrain the rest-frame K band magnitude of the AGN hosts, Mk^ for a comparison with 
the local M-qh-Mk scaling relations (Marconi and Hunt 2003; Graham 2007). Then, we proceed to 
derive the confidence interval on the measure of the total stellar mass. 

As a first step, we verified that, given our choice of AGN template, the derived K-band rest 
frame magnitudes of the host galaxies are not sensitive to the particular set of galaxy SED templates 
used (P07 vs. BC03). The average 1-a errors in the estimated Mk depend on the model SED used 
to fit the data, and are approximately 0.4 and 0.3 magnitudes for BC03 and P07, respectively. 
Keeping this in mind, we find a good agreement between these two methods: the difference in Mk 
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between the estimates obtained using the P07 templates and those obtained using the BC03 ones 
is strongly peaked at around zero, with small scatter (of the order of 0.3 dex), with the exception 
of 8/89 outliers (i.e. objects with |AM/f | > 0.7, larger than the sum of typical 1-a errors in Mk)- 
The estimated Mk (in VegJ^. for the fits with the BC03 templates) are given in Table [H while 
the distribution of the estimated Mj^ is shown in the upper panel of Fig. [2l There, as a term of 
reference, we have also plotted the corresponding distribution of the rest-frame K-band magnitudes 
of the AGN components from our fits. Typically, AGN are 1-2 magnitudes brighter than their hosts 
in the K-band; thus, the fraction of K-band light contributed by the host galaxy has a median of 
about 25-30%, with 90% of the object having this ratio smaller than 0.4. The ratio /gai,K = 
of the galaxy-to-AGN luminosity in the rest frame K-band is also given in table [TJ 

Encouraged by the robustness of the determination for the AGN hosts in our sample, we 
proceed to the discussion of the stellar mass estimates, and refer the reader to section 16.31 for a 
further discussion of how a different choice of AGN SED template could modify our results. 



3.2. Host galaxy masses 

The fitting procedure with the BC03 templates allows us to estimate also the total stellar 
mass of the AGN hosts: each combination of SFH, r and tage (see § [3|) is uniquely associated to a 
value of SFR and total stellar mass, computed taking into account the effect of stellar mass loss. 
The distribution of these stellar mass measurements is shown in the right panel of Figure [2l The 
estimated LogM* (for the fits with the BC03 templates) are also given in Table [TJ 

What kind of galaxies are these? As we will discuss in more detail later on, it is extremely 
difficult to extract reliable information on the star-forming properties of these objects, due to the 
dominant presence of the AGN emission. A more statistical comparison can be made by comparing 
the inferred masses with the overall galaxy mass function in the same redshift range as obtained 
by the S-COSMOS survey (duly shifted in mass to account for the difference in the average stellar 
mass between the Chabrier and Salpeter IMF, see Ilbert et al. 2009). We do not detect any galaxy 
with mass larger than lO^^'^M©, and most of our objec ts have stellar mas ses between 10^'^'^ and 



IO^^-^Mq. According to the S-COSMOS mass functions (jllbert et al. \\200% . within this range, 
only less than ~ 40, 35 and 20 % of all galaxies are "quiescent" (i.e. lie on the red sequence) for 
redshift ranges [1, 1.2], [1.2, 1.5] and [1.5, 2.0], respectively. 

A similar conclusion about the average star-formation properties of the type-1 AGN hosts 
can be drawn by studying the mass to light ratio of the best fitting models. Figure [2] shows the 
relation between rest frame K band luminosity and total stellar mass. Once again, we can compare 
these estimates of the M/Lk ratio with those measured by Ilbert et al. (2009) in the COSMOS 



Since our magnitudes were all calibrated on the AB system, in order to ease the comparison with literature work, 
we use the following conversion between Vega and AB COSMOS K-band magnitudes: Mx.voga = Mk, ab — 1-84 
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IR-selected galaxy sample (see Ilbert et al. 2009, appendix D). Taking into account the f»0.24 
dex shift due to the different IMF choices, our mean Log(M/L/f) ^ —0.55 follows more closely 
the expectations for the so-called blue-cloud (star- forming) galaxies, but the nominal uncertainties 
remain largj^. 

As a further test of our method, we have compared the total stellar mass estimates with those 
derived by Jahnke et al. (2009) for a small sample (18 objects) of X-ray selected AGN in the 
COSMOS field for which simultaneous HST/ ACS and HST/mCMOS observations allow to derive 
mass to light ratios (and stellar masses) of the resolved hosts. We find that the two independent 
methods are in broad agreement with each other. 5/18 of them have also zCOSMOS spectra, and 
are part of the sample described here (see the empty stars in Fig. [5]). Their total stellar masses 
calculated with the two methods agree to better than 0.2 dex, with our estimates being consistently 
on the lower side. We discuss this issue further in section 17.11 



4. Virial black hole mass measurements 

Due to the uncertainties in the actual geometry of the broad line region of type 1 AGN (see e.g. 
Mc Lure and Dunlop 2001), and/or to the different choice of absolute calibration of the black hole 
masses of local AGN studied with reverberation mapping, there exists a large number of different 
formulae in the literature that relate black hole mass, Mgll line width and continuum luminosity 
(see e.g. Mc Gill et al. 2008, and references therein). They can all be expressed in the form: 

log ^ =A + log{FWHMfQ^Q{XLsoooA4f) (1) 

where FWHMiqqq is the FWHM of the line in units of 1000 km s~^ , and AL3000.44 is the contin- 
uum luminosity at 3000A in units of 10^^ erg s^^. For Mgll lines, Mc Lure and Dunlop (2004) 
proposed A=6.51 and (3 = 0.62, based on a fit to the radius-luminosity relation for AGN with 
reverberation mapping and total luminosity XLx > 10^^ erg/s, while Vestergaard et al. (2009; in 
preparation; see also Trump et al. 2009) propose a different scaling (A=6.86, (3 = 0.5). Each 
of these relations carries a significant scatter of about 0.3 dex (McGill et al. 2008). The expo - 



nent (3 is related to the empirical calibration of the radius-luminosity relation (jKaspi et al. II2OOOI ) 



A recent re-ana lysis of reverberat ion mapping observations, fully accounting for host-galaxy light 



contamination (jBentz et al. II2009I ). points towards a value of 0.45 < /3 < 0.59, consistent with the 
simple expectations from photoionization models of the broad line region (/3 = 0.5). Here we will 
adopt the relation derived by McGill et al. (2008), with A=6.77 and (3 = 0.47. This was derived 
by cross-calibrating a number of different estimators applied to BLAGN in a redshift range where 



^® We note here that the level of uncertainty of the SFR for our AGN hosts is such that we cannot clearly discriminate 
between what, in the galaxy formation jargon, is usually called "blue cloud" and "green valley" (see e.g. Silverman 
et al. 2008) 
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more than one broad emission line can be observed simultaneously in optical spectra. Additional 
systematic errors in the derived black hole masses introduced by the use of any of these relations 
will be further discussed in section (see also the discussion in Treu et al. 2007). 

We do not take into account possible effects due to radiation pressure on the Broad Line Region 
(BLR), that could lead to systematic under-estimate of the black hole mass for the objects with 
the higher Eddington ratio. Such an effect has been estimated empirically by Marconi et al. (2008) 
for the H/3 broad emission in a local AGN sample. No calibration is currently available for Mgll 
emission lines. Moreover, the size of the possible radiation pressure correction is expected to scale 
linearly with the inverse of the Eddington ratio, and should not be too large for the moderate- 
luminosity AGN in our sample (see section [^31 below) . We thus made the conservative choice not 
to include any radiation pressure correction to the virial relation ([T]). As we will discuss at length 
in the following, this corresponds to minimize the amount of possible evolution detected in the 
average M-qr/M^ ratio as a function of redshift. 

In the following subsections, we will describe the method used to measure both line width and 
continuum luminosity needed to apply eq. ([1]). The measured FWHMiqqq, L300O) together with 
their la errors, are given in Table [2j 

4.1. FWHM measurement 

Type— 1 (un-absorbed) AGN spectra in the wavelength region of interest are usually character- 
ized by a power-law continuum, of the form /a oc A", broad line emission from Mgll plus a complex 
of Fell emission lines, also broadened at the typical velocities of the broad line region (Boroson 
and Green 1992; Vestergaard &; Wilkes 2001). Accurate subtraction of the broad Fell features 
is thus an important step in the process of obtaining broad Mgll line widths. It is not always 
straightforward to keep track of the Fe II subtraction technique in previous studies high-redshift 
scaling relations. For example, Peng et al. (2006a) erroneously report that Mc Lure and Jarvis 
(2002) did not performed any Fell subtraction, contrary to what stated in section 3.2 of Mc Lure 
and Jarvis (2002). For the sake of completeness, we briefly discuss in section [6^2] what would be 
the typical systematic effect of not removing an Iron template from the spectra. 

We adopt here the theoretical Fell template calculated by Bruhweiler & Verner (2008), 
calibrated by fitting the Seyfert 1 galaxy IZwl spectrum (model Fe_dll-m20-20.5 available at 
http : //iacs . cua.edu/people/verner/FeIIf). In order to apply the iron template to the spectra, 
the line width of the template must be matched to that of the AGN spectrum. We achieve this by 
means of an iterative procedure, as described below. 

First of all, we convolved the original Fe II template with gaussian functions of different widths, 
ranging from v = 1000 km s~^ to v = 15000 km s~^ in step of 250 km s~^ to produce a grid of 
broadened templates. From the original AGN spectrum we then derived: (a) a first rough estimation 
of the FWHM of the Mg II line using a single gaussian fit with iraf-splot package and (b) the mean 
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values of the continuum flux and the corresponding errors (derived from the noise spectra) in the 
following wavelength windows: [2660A-2700A]; [2930A-2970A]; [2715A-2750A]; [2850A-2885A] and 



described by a combination of power law and Fe II emission smeared by a velocity width v equal 
to the one derived from the Mgll emission line 



In order to find the best fitting model, we then performed a chi-square minimization and we thus 
derive the three parameters a, a and h. Finally, we subtracted the best fit power-law plus Fell 
model to the original spectrum obtaining a new spectrum that contains only the Mgll emission 
line. This is done using the iraf-sarith package. 

We then measured the FWHM of the line, modeled using one, two or three gaussians (ab- 
sorption and/or emission) and we choose the best solution according to best reduced computed 
m the [2650A-2950A] range. The best fit model for the emission line spectrum is thus given by: 
Mx = X^^o 0'iGx{\eak,i, o'i), where N is the number of gaussian components chosen by the fit. For 
each gaussian component Gx, a-i is the intensity (positive for emission, negative for absorption), 
\eak,i the Wavelength of the peak and cjj the width. Each of these quantities carries a statistical 
error provided by the splot-iraf package. 

We then compute the FWHM of the model emission line complex, and the error on the total 
FWHM, 

^ fwhm is computed by propagating the errors on the single gaussian components of the fit. 

Using the new FWHM determination for the line, we then iterated the whole procedure (we subtract 

the new broadened iron template chosen according to the FWHM measured from the original 

spectrum), until the fits converge and the final measure of the FWHM is stable. The final FWHM 

measurement is corrected for the finite spectrograph resolution assuming that FW H M^^^^-^^-^ = 
x2 

FWHM^^^ - where R is the mean instrumental resolution that for the zCOSMOS spectra 
is ~580. The best fit decomposition of the Mg II region of the spectra is shown in Figure [3] for 
the same nine AGN whose SED decomposition is shown in Fig. [TJ For three objects the fitting 
routine requires the presence of absorption in the Mgll line region (one is shown in the top left 
panel in figure [3]) . They are marked as blue open circles in figure Some residual errors due to 
over-subtraction of the iron template red- wards of the Mgll emission are apparent in some cases. 
It is clear that a single Fell template may not be adequate to fully describe this complex emission 
in all objects. A detailed study of the properties of the Fell emission requires higher signal-to-noise 
spectra, and is beyond the scope of the present paper. 




|. We assumed that in these wavelength ranges, the spectrum can be completely 



F{\) = aA° + 6 X FeII{X,v). 



(2) 



•^"For a handful of objects (5) with relatively small FWHM, a significant improvement in the fit is obtained by 
shifting the two central windows for continuum plus Fell fitting to [2730A-2770A]; [2830A-2870A], closer to the most 
prominent emission peaks of the Fell complex. 
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4.2. 3000A and bolometric Luminosity 

From the best fitting power-law continuum we derive an estimate of the AGN luminosity at 
3000A. Our points have been corrected for aperture effects normalizing the i-band fluxes measured 
in the VIMOS spectra with the observed ACS i-band photometry. Monochromatic continuum 
luminosities at 3000A were then calculated from the average best fit continuum flux rescaled in the 
2980-3020 A rest frame. The error on the continuum luminosity is obtained from the average of 
the noise spectrum in the same wavelength range. 

From the measured total AL3000 AGN luminosity we derive the bolometric one using the 
luminosity-dependent bolometric correction factor /boi of Hopkins et al. (2007). The distribution 
of bolometric luminosities for the type-1 AGN in our sample is shown in the upper panel of Figured 
Assuming a standard radiative efficiency of erad = 0.1, the median accretion rate onto the black 
holes, M = -/^boi/eradC^) is of the order of 0.4 Mq/jt:. 



4.3. Black hole masses and Eddington ratio distribution 

We have tested the impact of the choice of single-epoch virial formula for the black hole mass 
uncertainty. The typical spread in LogMBH among the three different estimators discussed in sec- 
tion [3] is of the order of 0.2 dex, similar to the observed scatter in the virial relations themselves 
(Vestergaard & Peterson 2006; McGill et al. 2008). In the following, we fix the systematic un- 
certainty in the LogMen determination to 0.2 dex, and for our statistical analysis of the scaling 
relation evolution, we assign to each black hole mass measurement an error given by the sum of 
the statistical and systematic uncertainties. 

Figure H] shows the estimated black hole masses vs. the bolometric luminosities for all the 
objects in our sample (upper panel) as well as their location in the Luminosity-Eddington ratio 
(A = iboi/-^Edd) where L^dd = 1-3 x W^^Mbh/Mq is the Eddington lumino sity) plane in the lower 



panel. As typical for optically selected sa mples of type-1 AGN/QSOs (jKollmeier et al. 1 120061 : 



Gavignaud et al. 1120081 : iTrump et al. II2009I . and references therein), the distribution of Eddington 



ratios is quite narrowly distributed, with a median of ~ 0.1. 

Analogously to what found in BLAGN samples selected in similar redshift ranges and at 
comparable depths (VVDS, Gavignaud et al. 2008), we also observe a trend of increasing Eddington 
ratio with increasing bolometric luminosity of the AGN (see also Netzer et al. 2007, for a sample 
of higher-redshift QSOs). Fitting a straight line we find LogA oc 0.64LogLboi5 but the slope of 
such a relation depends critically on the exponent /3 in the adopted virial relation ([T|). As already 
pointed out in Gavignaud et al. (2008), if we chose the McLure and Dunlop (2004) formula (with 
/3=0.62), we would obtain a shallower slope: LogA oc 0.47LogLboi- Selection effects could certainly 
be playing a role in determining the distribution of sources in the A — Lboi plane (for example, 
objects in the lower right corner, i.e. massive black holes at low accretion rates, probe the massive 
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end of the SMBH mass function, which is rapidly decUning in this redshift range, see e.g. Merloni 
&; Heinz 2008). The detailed distribution of AGN lifetimes as a function of luminosity can also be 
responsible for the observed trends, as suggested by some numerical models for QSO lightcurves 
(Hopkins and Hernquist 2009). For further discussion of possible selection and/or systematics 
effects, we refer the reader to the work of Trump et al. (2009), who have studied in greater detail 
the larger sample of XMM-COSMOS AGN with IMACS spectroscopy. 

5. Scaling relations and their evolution 

In this section, we quantify the amount of evolution (if any) in the scaling relations between 
nuclear black holes and host galaxy properties observed locally. 

We begin by showing in the left panel of Figure [5] the location of our AGN (filled circles) in the 
Log M^ii-Mx plane. As a reference, we show there the best fit relation derived from local inactive 
galaxies by Graham (2007) and given by 

LogMBH = 8.29 - 0.37{Mk + 24) (3) 

(to guide the eye, we show also the same relation offset by ±0.33 dex, the total scatter in the 
Graham (2007) relation). Our objects are shifted towards brighter hosts and/or smaller black hole 
masses with respect to the z = relation. The objects do not seem to obey any tight relation 
between black hole mass and Mk, although the range of Mbh probed is relatively limited. If one 
tried, for the sake of comparison, to fit the sample with a relation with the same slope as in eq.([3]), 
the best fit normalization would be shifted by ~ 0.28 dex, and the intrinsic scatter would be as large 
as 0.44 dex. We have also measured the offset, Tk{z) (here defined as the projected distance in the 
LogA/BH-LogLx plane), of each observed point from the local relation and studied its evolution as 
a function of redshift. The best fit obtained imposing the functional form Tk{z) = 5iLog(l + z), 
gives 5i = —0.73 ±0.08 (these results are unchanged if we instead adopt the Mk derived by fitting 
the SED with the P07 template). The fit are performed taking into account both errors and lower 
limits on F, using a Monte Carlo approach, described in Bianchi et al. (2007). 

Such an evolution is somewhat stronger than that observed by Peng et al. (2006b), where 
the R-band luminosity of 30 lensed and 20 non-lensed QSO hosts was measured based on detailed 
HST image modeling of 1 < z < 4.5 quasars. Their QSOs apparently lie almost exactly on the 
same observed Mbh--^^r relation as their z = relic counterparts. However, as already noted by 
Peng et al. (2006b), high redshift galaxies shall have a different mass to light ratio as compared 
to their z = descendants, at the very least because of passive evolution of the stellar population, 
in the extreme case of non-star forming galaxies. This can be simply accounted for by passively 
evolving the rest-frame K-band luminosities down to z = to allow a more direct comparison with 
the local relation. We have done this using once again the Bruzual and Chariot (2003) template 
libraries, and assuming all objects formed all their stars in burst at zj = 3. Their local descendants 
(passively evolved) would have dimmed to the magnitudes indicated by the black empty circles in 
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the right panels of Fig. [5l which are now broadly consistent with the local Mbh — Mk relation, 
but slightly offset towards large black hole to host galaxy ratios. Such a shift towards a positive 
offset from the local relation would in fact be even stronger if indeed the host galaxies of our AGN 
sample were dominated, at least statistically, by actively star forming galaxies, as discussed earlier 
in section I 



This is indeed what happens when we consider the relationship between measured black hole 
masses and host total stellar masses, , obtained from the BC03 fits. In the right panels of FigureO 
we show the location of 89 zCOSMOS AGN in the LogMeH-LogM^, plane. In these panels, we have 
separated the objects on the basis of the measured /gai,K host-to- AGN luminosity ratio in the rest 
frame K-band (with open triangles having /gai,K < 0.34 and filled squares /gai,K > 0.34). As a 
ref erence, a solid hne sh ows the local best fit relation between black holes and spheroids as derived 



bv lHaring fc Rix I tOM ): 

LogMBH = -4.12 + 1.12(LogM,) (4) 

Our objects now show a modest, but clear, offset from the local relation, especially at the high 
black hole mass end (see also Figure [8]) and in the highest redshift bin (upper right panel of Fig. [5]). 
Radio-detected AGN (squares) appear to be distributed similarly to the rest of the population, but 
low number statistics prevent us from reaching any firmer conclusion on their properties. 

As our method is affected by substantial uncertainties in both black hole and host galaxy mass, 
and none of the two can be treated as truly independent variable, we chose to measure the deviation 
of our data-set from the local scalin g relation by measuring the distance ALog(MRH/M*) of each 



point from the Haring &: Rix (2004) relation, perpendicular to the relation itsel^J- 



Figure El shows the measured offset of each point from the local relation as a function of redshift. 
The black solid line shows the best fit obtained imposing the functional form ALog(AfBH/-^^*)(-2^) = 
52Log(l+2:), where we find 62 = 0.68=b0.12. Also in this case we have used a Monte Carlo simulation 
to derive the best fit linear regression coefficients; lower limits were treated as if the true value of 
ALog(MBH/Af*) were uniformly distributed up to a common value of 1.2. 

The inset of Fig. [6] shows our dataset and best fit evolution, in black, as compared to a 
number of estimates at lower or comparable redshift. The best-fit evolution from the zCOSMOS 
data is in reasonable agreement with previous estimates both at lower (Salviander et al. 2007) 
and at higher redshift (Peng et al. 2006b). We note here that the significant amount of scatter 
in our dataset translates into a relatively weak statistical significance of the measured offset (see 
the binned points in the inset of Fig. [6]). Moreover, any redshift dependence within our sample 
only is not statistically significant. Although weaker, the evolution we measure is also marginally 



■^^Given the slope in the lHaring fc Rix I (2004) relation A = 1.12, the measured offset multiplied by S = y/1 + ~ 
1.5 gives the increase in black hole mass AMbh given a host galaxy mass, compared to the local value. This should 
be kept in mind when comparing with results from previous works, which usually measure the offset from the scaling 
relations in terms of "excess black hole mass", i.e. vertically in Fig. [5] 
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consistent with that observed by Treu et ah (2007) and Woo et al. (2008) from their sample of 
Seyfert galaxies at z = 0.36 and z = 0.57. 

It is interesting to notice here that the majority of observational data-points for AGN sample 
do indeed show a broadly consistent amount of offset at all redshifts probed. This might suggest 
that intrinsic differences in the SMBH/host galaxy relation between active and inactive galaxies 
could play an important role besides any genuine cosmological evolution. Large, uniformly selected 
samples of AGN hosts, spanning a larger redshift range than probed here (e.g. Decarli et al. 2009), 
as well as accurate comparisons of scaling relations for active (reverberation mapped) and non-active 
galaxies at low-z (see e.g. Onken et al. 2004), are and will be very important in disentangling true 
redshift evolution from other systematic differences with the local samples. 

Finally, no significant trend is found by dividing our sample into fast and slow accretors on 
the basis of their measured Eddington ratio, contrary to the results of a number of studies of local 
AGN, which have found evidence that high-accretion rate objects have smaller Mbh/-^-^* ratio as 
compared to less active AGN (Greene and Ho 2006; Shen et al. 2008; Kim et al. 2008). 



5.1. Flow patterns in the Mbr-M^, plane 



As we have mentioned before, constraining the star formation rate of the AGN hosts in our 
sample is a very difficult task, given the strong AGN component dominating the SED in the rest- 
frame optical/UV bands. Individual estimates of SFR based on the minimization procedure 
described above give results which are uncertain by up to 0.7-0.8 dex. Nevertheless, despite these 
very large uncertainties on the estimated star formation rates for our AGN hosts, we can try to 
assess the general direction of motion of the objects in the Mbh — Mk plane; although for each 
individual object it will be hard to accurately pin down the change in total stellar mass, the 
overall "ensemble" average motion of the flow could give interesting indications on the longer term 
evolution of the scaling relations. 

To this end, we show in Fig. [7| as red thin arrows the predicted flow patterns of the BLAGN. 
The tip of each arrow marks the location where the system will find itself within tstar = 300 Myr, if 
continually forming stars at the estimated star formation rates, while at the same time the central 
black hole keeps accreting at the measured rates for a fraction of this time equal to the AGN 
duty cycle 6t. The exact value of igtar is of course arbitrary, and has been chosen in order to 
better visualize the flow pattern. Changing it, will simply rescale the length of the arrows, but 
won't change their orientation. For each object observed at redshift z having nuclear bolometric 
luminosity Lboi and host galaxy stellar mass M^,, the duty cycle is estimated by ta king the ratio 



of the AGN bolometric luminosity function (IHopkins. Ri chare 



the mass function of highly star-forming galaxies (jllbert et al 



s fc Hernquist 1 120071 ) (/>(Lboi,-z) to 



20091 ). which are assumed here to 
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represent the parent population of the sample at hand, i;^gai,HSF(-^-^*) -z)- 

5t(Lboi,M,,z) = — — - (5) 

(Pga.l,HSF [Jyh,Z) 

As the mass functions have a sharp exponential cut-off at high masses, as opposed to the power-law 
decline of the QSO luminosity function, the duty cycle defined above increase quickly for the most 
massive hosts, causing the upwards turn of the flow pattern. Intriguingly, the predicted motion of 
the objects in the Mbr-M^ plane does lead to a reduced scatter of the points, i.e. the flow appears 
to be "converging". To give a simple quantitative estimate of this effect, we measure the scatter 
by fitting all points (apart from those with only upper limits for the host galaxy properties) with 
a linear relation. Given the reduced dynamic range in black hole masses we are probing, and for 
the sake of simplicity, we fix the slope of the correlation to the locally measured one (1.12; Haring 
and Rix 2004), and let the normalization be a free parameter. The observed points have thus a 
normalization of about —3.7 and an intrinsic scatter of 0.43. The tips of the arrows, instead, move 
closer to the local relation (normalization ~ —3.9) and with a much reduced scatter of 0.34. 

Future observations of lower redshift AGN, for which an accurate determination of both to- 
tal stellar mass and star formation rate is more easily achieved by a combination of multi-band 
photometry and high-resolution spectroscopy will provide much better (and more reliable) maps of 
the flow patterns of the AGN-host galaxy systems, revealing fundamental details on their physical 
coupling. 



6. Systematics and selection effects 



The main result of the previous section is that our estimates of the type-1 AGN host physical 
parameters are (although marginally) inconsistent with the hypothesis that they lie on the z = 
scaling relation. Here we wish to discuss how much of this observed offset can be due to various 
systematics and selection effects. We identify here two kind of biases: the first is the combination 
of systematics inherent to our methods to measure either BH or host galaxy's mass, and we will 
discuss these first (sections 16.11 16.21 and 16. 3p . The second, more subtle, is a 'luminosity /mass func- 
tion weig hted' bias on any AGN-selected sample introduced by the intrinsic scatter in the scaling 
relations dAdelberger fc Steidel II2OO5I : IPine et a:. Ibood : iLauer et al.lboOTl : iTreu et al. Iboo?! ) . This 
could induce a spurious effect on the measured offsets, provided the scatter in the relation is large 
enough. We will discuss this effect in some detail in section | 



6.1. IMF and galaxy stellar masses 

In this work, mainly to allow a direct comparison with most previous works on the subject, 
we have adopted a Salpeter IMF to calculate the total stellar mass of the host galaxies based on 
the BC03 SED fitting procedure. As mentioned already in section [3l adopting a different IMF 
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will result in a systematic shift of the best values for M,,,, typically reducing the stellar mass, and 
increasing the ratio Mbh/A^*. To quantify the systematic uncertainty in the measured evolution 
introduced by the uncertainty in the IMF, we have re-calculated the stellar masses using a Chabrier 
IMF. We thus shifted the estimated values of the total stellar mass of the AGN hosts by -0.255 
dex (Pozzetti et al. 2007) and found the following values for the exponent of the redshift evolution 
function: 82 = 1-15 it 0.13. The smaller host masses implied by the new choice of IMF result in 
a larger positive offset of the Mbh/-^* ratio from the locally determined value, requiring a more 
pronounced evolution. On the other hand, it should be kept in mind that realistic SFH can be 
different from the smooth ones adopted here, due to the presence of short bursts of star formation. 
In fact, Pozzetti et al. (2007) have tested how the measured stellar masses changes when random 
bursts are superimposed on smooth SFH (with a similar range of combinations of r and tagc to the 
one chosen here). They conclude that, for z > 1.2 samples, the mean mass computed with smooth 
SFH is on average 0.16 dex smaller than the one computed using a bursty SFH (see fig. 5 in Pozzetti 
et al. 2007). This offset is smaller than our typical errors on the stellar masses. Moreover, a robust 
assessment of the "burstiness" of the SFH cannot be performed on our sample, as the number of 
extra parameters required (duration of a burst, fraction of stellar mass produced in the burst, time 
since last one) would surely introduce strong degeneracies in our fits. 

Another well known source of systematic uncertainty lies in the choice of the BC03 templates 
to describe the stellar populations of the AGN hosts. Different groups have in recent years built 
different stellar population models using different treatments of stellar structure and evolution (see 
e.g. Fioc &: Rocca-Volmerange (1997), Silva et al. (1998), Maraston (2005)). A thorough discussion 
of the relative differences among them is far beyond the scope and the aims of this paper (but see, 
e.g., Longhetti and Saracco 2009). Suffice it to say, in this context, that the typical offset in the 
estimated masses of test galaxies with ages similar to those considered here are much smaller (of 
the order of ±0.1 — 0.15 dex, see e.g. Cimatti et al. 2008; Longhetti and Saracco 2009) than the 
statistical uncertainties of our mass measures. 

Finally, it is clear that by comparing total stellar masses of AGN hosts with the bulge/spheroid 
masses of the local galaxies originally used to derive and identify the scaling relations, we are 
introducing a significant bias (see also Jahnke et al. 2009; Bennert et al. 2009). Lacking any 
imaging information and reliable bulge-to-disc (B/T) decomposition, we can just argue that, at 
the very least, the black hole to bulge mass ratio should show an even larger offset from the local 
scaling relation. Systematic trends in the B/T ratio with redshift (see e.g. Merloni, Rudnick and 
Di Matteo 2004) will need also to be taken into account. 



6.2. Black hole mass measurements 

We have already mentioned in section [J] that there is currently a substantial uncertainty on 
the actual parameter of the virial relationship to be used for the estimate of black hole masses in 
broad line AGN. In the calculations so far, we have adopted the McGill et al. (2008) expression, 
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but different relationships, based on different calibrations and/or assumptions about the broad hne 
region geometry exist in the hterature. 

To test the systematic effects on the measured evolution of the black hole to host galaxy mass 
ratio evolution, we have re-calculated black hole masses using the MLD04 and V09 relations (see 

§[D- 

Adopting the relation of McLure and Dunlop (2004), which has a steeper BLR size-luminosity 
relation (0.62 instead of 0.5), substantially reduces the amount of observed evolution: we obtain, 
for the exponents of the redshift evolution function, 82 = 0.47 it 0.12, bringing it closer to the 
expectations of a purely luminosity-bias dominated effect if the Mbh — relation has an intrinsic 
scatter as large as 0.5 dex (see section below). On the other hand, adopting the V09 relation, which 
produce, on average, larger black hole masses than in our fiducial case, we obtain 62 = 0.91 ± 0.12, 
indicating a larger amount of positive evolution. It is worth keeping in mind, however, that the 
McLure and Dunlop (2004) normalization of the "virial" black hole mass estimate is based on a 
specific (theoretically motivated) assumption on the BLR geometry, while the V09 one is empirically 
calibrated with the local Mbh — relation, i.e. it does not assume any preferred geometry. 

Other possible systematic effects could be due to evolution of the physical properties of the 
broad line region itself. However, no significant trend with redshift was found in our data for 
either the Fell emission strength (measured relative to the continuum luminosity), or the overall 
goodness of the continuum fits, that could have signalled an inadequacy in our model of the Fell 
template (due, for example, to systematic changes in metallicity) . Higher signal-to-noise spectra 
are probably needed to assess these issues with the due care. 

Finally, we would like to note here that, when comparing with results present in literature, 
one should take into account the different methods of line fitting and Fe II subtraction. A thorough 
analysis of the induced bias is not straightforward, as the specific techniques of broad line fitting 
can be different. Nevertheless, we have ourselves performed a new fit of the Mgll line complex 
without including any Fe II template in the continuum. The distribution of the (log of the) ratio 
between the FWHM calculated without and with Fell subtraction is centred at a positive offset 
of about 0.1 dex (which would correspond to a difference in mass of about 0.2 dex), but with a 
significant tail of higher ratios towards low values of the FWHM. 

6.3. Choice of AGN SED 

In sectionOwe have shown how the measure of the rest-frame K-band luminosity of the BLAGN 
hosts does not depend strongly on the set of galaxy SED templates used to fit the composite 
AGN+galaxy spectral energy distribution, provided that the same AGN SED is used (we adopt 
here the Richards et al. (2006) mean QSO SED, allowing for additional dust extinction of the 
nuclear light). 
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An obvious possible objection to our results is that, if the AGN SED were markedly different 
from the Richards et al. (2006) template at the typical luminosities of the zCOMSOS bright sample, 
we would be introducing a severe bias in our measures of the host galaxy parameters. More 
worryingly, any systematic trend whereby AGN SEDs change with either redshift, luminosity, 
black hole mass or Eddington ratio, would introduce spurious trends in the measured evolution of 
the offset from the local scaling relations. 

Elvis et al. (2009, hereafter E09) performed a thorough, systematic study of the spectral energy 
distribution of X-ray (XMM-COSMOS) selected, spectroscopically confirmed, type~l (broad line) 
AGN in the COSMOS field. The interested reader is referred to E09 for a discussion of the main 
properties of the sample. Here we would like to point out that, even when selecting only those 
objects that are classified as "pointlike" from their ACS images, in the redshift range of interest 
here (1 < z < 2.2) a non-negligible contribution due to the compact stellar emission from the host 
is still significantly present in the observed SED, even in the brightest sources. 

Indeed, the mean SED of the "pointhke" type-1 COSMOS AGN with Lboi > lO^^'^ shows 
a less pronounced dip in the 1/im region than either the Elvis et al. (1994) or the Richards et 
al. (2006) templates. As a test, we have used this new, COSMOS based, AGN SED together with 
the BC03 templates to fit the composite spectral energy distribution of the objects in our sample. 
As expected, the flatter mean AGN SED results in a lower residual host galaxy contribution. 
The number of "undetected" hosts (i.e. they have only upper limits in the rest frame K-band 
magnitude estimate) rises from 10 to 24; on average, the hosts result about 0.4-0.5 magnitudes 
fainter. Consequently, the positive offset of type-1 AGN from the local scaling relations would be 
larger (we have measured, for these choice of AGN SED, 82 — 1.3). 

We argue that, on the basis of the photometric data only, it is not possible to decide whether 
the apparent flatness of the AGN SED at NIR wavelengths is indeed due to a dramatic change 
of intrinsic nuclear continuum emission (thus minimizing the host galaxy contribution), or, on 
the contrary, it is mainly due to stellar light contamination. As the QSO SED we have adopted 
represent an extreme within the COSMOS database (see E09), our approach, which maximizes 
the host's contribution is the most conservative one with respect to the measured evolution of the 
Mbh/-^* ratio. Independent clues on the intrinsic nuclear continuum shape from e.g. detailed 
spectroscopic and/or polarimetric studies will be helpful. 

Finally, we have als o considered th e effec ts of AGN variability on the observed SED. This has 



been studied in detail by lSalvato et al. I (120091 ). where it was demonstrated that correcting for even 
relatively small variation of the photometric points due to the non-simultaneous observation times 
can signiflcantly improve the photometric redshift determination for AGN. A "variability corrected" 
photometric catalogue is available for the subsample of X-ray detected type-1 AGN (82/89 objects 
in our sample). We have thus recomputed the hosts total stellar masses and studied the evolution 
of the Mbh/-^* ratio also for them, flnding no significant difference in the redshift dependence of 
the measured offset (^2 = 0.84 ± 0.13). 



- 23 - 



6.4. 



Selection bias 



In this section we estimate the possible bias due to selection effects. Our objects are selected 
essentially on the basis of the nuclear ( AGN) luminosity, and on the detectability of the broad Mg II 
emission line, clearly leading to a bias towards more massive black holes, similar to Malmquist 
(1924) bias for luminosity selected samples of standard candles. 

In fact, there is a more subtle effect, generally applicable to all cases where two properties 
of a class of objects are known to be correlated with certain intrinsic dispersion, and one wishes 
to determine the probability distribution of one of the two quantities, having selected objects on 
the basis of measurement of the other one. For purely flux limited samples, this bias was already 
discussed and calculated by Kellerman (1964, see Appendix), for the case in which the two quantities 
were the spectral indices of AGN at two different radio frequencies. This was then generalized by 
Francis (1993) to the case of AGN spectral slopes in any two given independent bands. For the 
specific case of black hole and host galaxy masses (or velocity dispersion), this bias was discussed 
already in Adelberger and Steidel (2005), but has been scrutinized in depth in Lauer et al. (2007), 
which also recovered the Kellerman (1964) results for the specific case of a flux limited sample, and 
we refer the curious reader to the Lauer et al. (2007) paper for a more thorough discussion. 

Here we make use of the main analytic results of Lauer et al. (2007) and apply them to our 
particular selection criteria. The null hypothesis we put under test is that the local scaling relation 
between black hole and host galaxy mass, assumed here to be given by eq. (j4]), with an intrinsic 
scatter a^, does not change with redshift, neither in normalization and slope, nor in scatter. 

It is in fact the intrinsic scatter in the local relation, together with the observed shape of the 
mass and/or luminosity function of the selected objects in the appropriate redshift ranges that 
determine the bias. In a nutshell, cosmic scatter in the Mbh — relation implies that there is 
a range of masses LogM* it for each object of a given black hole mass Mbhj where we have 
assumed, for simplicity, a symmetric scatter in the relation. If the number density of galaxies is 
falling off rapidly in the interval LogM,, it a^, it will then be more likely to find one of the more 
numerous small mass galaxies associated with the given black hole, and therefore a larger ratio 
M-QYi/Mtf. Thus, given a distribution of galaxy masses (mass function (/)(LogAf*)), and provided 
that the scatter is not too large, the logarithmic offset of each point from the correlation (jj]), 
assumed to be held fixed to the local determination, is given by [cfr. Lauer et al. (2007), eq.(14)]: 



where 1/5 = 0.67 (see footnote EI] in section [5]), /i = LogMen and {A,B) = (1.12,-4.12) are 
slope and intercept of the relation We estimate the logarithmic derivative of the mass function 
dLogM using the mass function determination from the S-COSMOS galaxy survey (Ilbert et al. 
2009; masses have been recalculated to adjust to our choice of Salpeter IMF) in the same redshift 
range probed by our BLAGN. Figure [8] shows our objects in the offset-black hole mass plane, where 
solid lines mark the expected bias from eq. ([6]) in three different redshift ranges, in the case of 




(6) 
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= 0.3. Dashed and dot-dashed hnes, instead, correspond to the cases of = 0.5 and 0.7, 
respectively. This plot clearly shows that the offset we measure is in excess to what expected in 
the most extreme case of large intrinsic scatter in the local relation, estimated by Novak et al. 
(2006) to be less than 0.5 dex (see also Giitelkin et al. 2009). This is mainly due to the fact that 
the AGN black hole masses we measure at the depth of the zCOSMOS selection function are not 
extremely large, and correspond, according to eq. (j4]) to a range of host galaxy masses where the 
mass function is not falling off too steeply. On the other hand, our observations could be explained 
in terms of luminosity bias only if the scatter in the Mbh — M^, relation were as large as 0.5-0.7 at 
z>l. 

Yet another test is possible, however. When studying the redshift evolution of ALog(MBH/Af*), 
fitting its redshift dependence with a functional form (fj5]), one does effectively take the average 
of the offset over a range of AGN luminosities above a well determined selection cut in any given 
redshift bin. In this case, it can be shown (Lauer et al. 2007, eq. 25) that the average offset at any 
given redshift is given by: 

{ALog(MB„/A/.))(.) = • 

where ^{L,z) is the type 1 AGN luminosity function at redshift z and Lmin and Lmax sue the 
minimum and maximum luminosity of the AGN that can enter our sample at the same redshift, 
given our survey selection function. We have calculated this bias adopting the VVDS type 1 AGN 
luminosity function of Bongiorno et al. (2007), and the results are plotted as red lines in Figure [H 
For our type-1 AGN sample, which extends well below the knee of the type 1 luminosity function 
in this redshift range, the expected bias in the offset is almost constant with redshift, amounting to 
w 0.3 (0.1) dex for = 0.5 (0.3). This is consistent wit h the estimat e for A LogMRH of Lauer et 



al. (2007) based on the local AGN luminosity function of lBovle et al. I (|2000l ). and slightly smaller 
than the measured offset. On the other hand, a larger intrinsic scatter in the scaling relation at 
the redshift of interest (up to = 0.7, dotted lines in Fig. [6]) could indeed be the cause of the 
measured offset in our sample. 

A note of caution is however in place. Our measured offset is based on a bl ack hole mass 



"viria l" estimator that has been implicitly calibrated on the local Mbh — relation (jOnken et al 



20041 ). This is usually done with the justification that the true geometry of the broad line region is 
not known, which enters as a multiplicative factor in the virial relation. However, the reverberation- 
mapped AGN used by Onken et al. (2004) to carry out this normalization, could themselves be 
affected by the 'luminosity /mass function weighted bias'. If that were the case, and if the local AGN 
samples spanned a similar mass and luminosity range as ours, the bias would be much reduced, 
and possibly canceled out completely (as suggested also by the results of Kim et al. 2008). The 
selection function of the Onken et al. (2004) AGN is, however, far from being well understood, and 
a quantitative estimate of the true expected residual bias is beyond the scope of this work. As we 
discussed above, large, uniformly selected samples of AGN hosts, spanning a larger redshift range 
than probed here, as well as accurate comparisons of scaling relations for active (reverberation 
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mapped) and non- active galaxies at low-z with similar selection functions, are and will be very 
important in disentangling true redshift evolution from other biases and systematic differences 
with the local samples. 

Finally, we would like to mention the further aspect of what, in general terms, can be defined 
luminosity bias of the AGN, namely the fact that faint broad line AGN cannot be detected in bright 
galaxies. This is clearly an important issue to consider in our case, given this zCOSMOS sample 
has a larger fraction of faint AGNs. 

In order to assess this, we have divided sources in Fig. [5] and Fig. [8] according to the measured 
luminosity ratio in the rest-frame K-band, with open triangles (filled squares) marking the most 
(least) AGN-dominated sources. It is interesting here to point out that, indeed, high contrast (AGN- 
dominated) objects do seem to slightly bias the result toward high Mbh/-^* ratios. However, 
removing all such objects, the main results presented here of a redshift evolution of the scaling 
relation, are not only confirmed, but do appear to be strengthened, as AGN-dominated objects 
(open triangles) are the largest outliers also in the lowest redshift bin. 

7. Discussion and conclusions 
7.1. Our results 

We have used an AGN-l-host galaxy SED decomposition technique to infer the physical prop- 
erties of the hosts of 89 (moderately luminous: i.e. mostly sub-L*) type-1 AGN in the zCOSMOS 
survey. Thanks to the deep, intensive multi-wavelength coverage of the COSMOS field, the ob- 
served spectral energy distributions are sampled to such a degree that the decomposition technique 
works reasonably well. We are thus able to derive rest frame K-band magnitudes and total stellar 
masses for the majority of our sample (80-90%, depending on the choice of AGN SED). Noticeably, 
our method allows us to properly quantify the uncertainty in each measurement. 

The bulk of the sample of BLAGN hosts we have studied have total stellar masses in the range 
lOW.5_iQii.3j^^_ Both the derived mass-to-light ratios and the sample average star formation rates 
seem to suggest that they are moderately-to-highly star forming objects, in good agreement with 
the host properties of both type 1 and type 2 (obscured) X-ray and/or IR selected AGN, either in 
a lower (Kauffmann et al. 2003; Jahnke et al. 2004; Hickox et al. 2009; Silverman et al. 2009) and 
in a similar (Brusa et al. 2009; Ross et al. 2009) redshift range. Reassuringly, also the fits with the 
POT phenomenological templates indicate an overall preference for star-forming hosts over passive 
(elliptical) ones. Obviously, a more accurate determination of the individual star formation rates 
is hampered by the dominant contribution of the AGN in the rest-frame UV/Optical part of the 
spectra. 

The black hole masses derived with the the "virial" or "empirically calibrated photo-ionization" 
method are broadly distributed around LogMeH ~ 8.5. Interestingly, the mass range probed by 
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our sample corresponds to that where the most stringent constraints on the z = scahng relations 
are available. The Eddington ratios A of our objects have a mean of LogA ~ —1. We observe a 
clear trend between Eddington ratio and bolometric lumin osity, that could be indicativ e of some 



specific luminosity-dependent AGN lifetimes distribution ([Hopkins Sz Hernquist I l2009l ). but the 



underlying effect of various selection biases in determining this trend needs to be further assessed 



(|Trump et al. 1 120091 ) . 



For a large number of objects in the sample (68/89; but the number falls to 24/89 if we consider 
the errors) the measured black hole to stellar host mass ratio is positively offset from that predicted 
by the local Haring & Rix (2004) local scaling relation. Assuming a redshift-dependent evolution 
in the ALog(MBH/Af*) of the form (^2Log(l + z), we measure 52 = 0.68 it 0.12^q 3, where the large 
asymmetric systematic errors stem from the uncertainties in the hosts' IMF, in the calibration of 
the virial relation used to estimate BH masses and in the mean QSO SED to be used. 

The scatter in the measured offset is substantial at all redshifts probed, such that we could 
still be consistent with a lack of evolution in the scaling relation at the 2-a level, as shown by the 
inset of figure O There, it is also apparent that the majority of observational data-points for AGN 
samples in the published literature do indeed show a broadly consistent amount of offset at all 
redshifts probed. This might suggest that intrinsic differences in the SMBH/host galaxy relation 
between active and inactive galaxies could play an important role besides any genuine cosmological 
evolution. 

Jahnke et al. (2009) have independently computed total stellar masses for a small sample 
(18 objects) of X-ray selected AGN in the COSMOS field for which simultaneous HST/ ACS and 
HST/NICMOS observations allow to clearly image the resolved hosts. The two independent meth- 
ods for measuring stellar masses are in broad agreement with each other, with a dispersion well 
within the (large) uncertainties that characterize each of these methods. 5/18 of them have also 
zCOSMOS spectra, and are part of the sample described here (empty stars in Fig. [5]). However, 
the objects in the Jahnke et al. (2009) sample do not show any significant offset from the local 
scaling relation. This might be due to a statistical fluctuation (more so given that they tend to lie 
in the lower redshift range probed by our sample, where the offset from the local scaling relation 
is smaller), or may indicate a more serious issue with the different estimates of the mass-to- light 
ratio of the AGN hosts. Larger samples of HST-imaged AGN hosts at longer wavelengths (with the 
newly installed WFC3) will be extremely important to settle this issue and, in general, to improve 
the calibration of our SED-based method to estimated host galaxy masses for AGN. 

We have taken particular care in examining the effects of the bias inevitably introduced by any 
intrinsic scatter in the BH-host mass scaling relation into any AGN-selected sample, as ours. We 
conclude that our data cannot possibly be explained if type-1 AGN and their hosts at 1 < z < 2 
lie on a scaling relation which has the same slope, normalization and scatter as the locally observed 
one. On the other hand positive evolution of the average Mbh/-^* ratio (larger black holes at early 
times in unobscured AGN) or of the intrinsic scatter (or a combination of the two) are needed to 
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explain our results. 



7.2. Implications for theoretical models 



What are the implications of these findings for our understanding of the cosmological co- 
evolution of black holes and galaxies? Let us briefly discuss recent theoretical investigations on this 
issue and the corresponding predictions for the evolution of the scaling relations. 

One of the earliest semi-analytic models to incorporate the evolution of supermassive black 
holes and the associated feedback effect were those by Granato et al. (2001,2004). There, triggering 
of AGN activity is not directly linked to merger activity, but rather generically to the process of 
bulge/spheroid formation. The rate of star formation and black hole accretion is regulated by the 
starlight radiation drag, and consequently, in a typical system the ratio Mbh/-^^* is initially small 
and rapidly grows until the AGN feedback sweeps the remaining gas. QSOs and, in general, type-1 
AGN are thus associated to the final stage of bulge formation, and it is very hard to produce any 
positive offset from the local relation like the one we measure. 

This is however a problem common to all feedback models in which the black hole energy injec- 
tion is v ery fast (explosive). Ind eed, the first published predictions of merger-induced AGN activity 
models (IRobertson et al. II2006I ) indicated that, if strong QSO feedback is respon s ible for rapidly 



terminating star formation in the bulge (jPi Matteo. Springel fc Hernquist 



2005 



Springel et al. 



20051 ). then very little evolution, as well as very little scatter , is expected for the scaling relations. 



Howe ver, later works within the same theoretical framework (iHopkins. et al 



2007 



Hopkins et al. 



20091 ) have analyzed in greater depths the role of dissipation in major mergers at different redshift. 
Under the assumption that black hole and spheroids obey a universal "black hole fundamental 
plane" (BHFP), where Mbh oc M^a'^, they show how, in gas-richer environments (at higher red- 
shift), dissipation effects may deepen the potential well around the black hole, allowing it to grow 
above the z = Mbh — relation, to a degree marginally consistent with our results. However, 
in the same physical framework, the Mbh — o'* relation is almost independen t on re d shift, which 



would contradict the observational results of IWoo et al. I (120061 ) ; iTreu et al. I (120071 ) ; IWoo et al. 



(12003). 



Another, related effect was discussed in Croton (2006). There it was assumed that major 
mergers can trigger both star formation in a bulge as well as black hole growth, in a fixed propor- 
tion. However, bulges can also acquire mass by disrupting stellar discs, a channel that should not 
contribute to black hole growth. The relative importance of these two paths of bulge formation 
may lead to lighter bulges for a given black hole mass at high redshift, as disks have a smaller 
stellar fraction. A subsequent study of this and other dynamical process of disk-to-bulge trans- 
formation was included in the work by Fontanot et al. (2006) and Malbon et al. (2007). They 
also confirmed qualitatively the predictions of Croton (2006), but found a much smaller effect, at 
most a factor ~ 2 at z = 2 in the Malbon et al. (2007) work, and preferentially for small mass 
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black holes Mbh ^ 10^ Mq. Even more comp lex semi-analytic mod els including various flavors of 



AGN-driven winds and their feedback effects ([Fontanot et al. II2006I ) can lead to various degrees of 



positive redshift evolution of the average Mbh/-^* ratio (Lamastra et al. 2009). 

As a general rule, we observe that, following increasing observational efforts to study the 
evolution of scaling relations, semi-analytic models have become more sophisticated over the years. 
This increase of sophistication's has allowed more complex behaviours of the coupled black holes- 
galaxy systems over cosmological times. As it is expected, more complex models also lead to an 
increase in the predicted scatter, even though a clear theoretical study on the redshift evolution of 
such scatter is still missing (but see the recent attempts by Lamastra et al. 2009 and Somerville 
2009). 



Hints from hydrodynamical simulations, both of isolated mergers (I Johansson. Naab &: Burkert 



2009 ) and of relatively small cosmological boxes ( Goldberg Sz di Matteo 20081 ) do indeed show a 



large scatter in the instantaneous ratio between black hole accretion and star formation rates, 
similar to what was found here (see also Silverman et al. 2009), thus suggesting that on the 
relatively short timescales over which un-absorbed AGN/QSOs are visible the physical connection 
between black holes growth and galaxy formation must be complex, too. How this would impact 
on the statistical and evolutionary properties of the galaxy population as a whole, however, is far 
from clear. 

The results we have presented in ^ coupled with analysis of selection biases of ^6.41 would 
suggest that a greater effort should be made by theoretical modellers to include a more accurate 
and realistic study of the evolution of the intrinsic scatter in any scaling relations, as well as that 
of slope and normalization. 

The "increased scatter" hypothesis, as an explanation of the observed offset, and its inevitable 
consequence that a 'luminosity function weighted' bias plays a significant role in the observed 
evolution of scaling relations, could be stren gthened if recent cla ims of under-massive black holes in 



IR selected galaxy samples were confirmed (jShapiro et al. II2009I ). Indeed, we should expect that in 



samples selected purely on the basis of the host galaxy stellar mass, rather than on AGN properties, 
the intrinsic scatter in the scaling relation should produce a bias going in the opposite direction 
as those discussed above, depending on the exact shape of the BH mass function at the redshift 
considered. 



7.3. Concluding remarks 

By taking advantage of the unique combination of VLT spectroscopy and deep multi- wavelength 
coverage of the GOSMOS field, we have presented here a novel method to study the physical link 
between supermassive black holes and their host galaxies in type-1 (un-obscured) AGN in the 
crucial redshift range 1 < z < 2. The main focus of this work is on the capability of our SED 
decomposition technique to provide reliable estimates of the total stellar mass of the AGN hosts. 
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and, even more importantly, reliable estimates of its uncertainty. 

The main result of our study is the observation of an offset in the Mbh — relation, such 
that, in the redshift range probed, for their given hosts black holes are on average 2-3 times larger 
than their counterparts in the nuclei of nearby inactive galaxies. A thorough analysis of all possible 
observational biases induced by intrinsic scatter in the scaling relations reinforces the conclusion 
that an evolution of the Mbh — relation must ensue for actively growing black holes at early 
times: either its overall normalization, or its intrinsic scatter (or both) must increase significantly 
with redshift. 

We close with two recommendations for future studies of the subject. From the observational 
point of view, it will be very important to explore methods to derive robust black hole mass 
estimates in high redshift samples of obscured AGN, that can be selected purely on the basis 
of their host galaxy properties. Broad emission lines at longer wavelengths, where the effect of 
obscuration are less severe, could be very useful in this respect. Also, a better understanding of 
the differences in the hosts' properties of active and inactive black holes is needed to allow a more 
meaningful comparison with the local scaling relations, and a better assessment of their evolution. 
From the theoretical point of view, more efforts should be devoted to derive robust predictions for 
the coupled evolution of slope, normalization and intrinsic scatter in the scaling relations, and to 
properly include in the models the selection effects that clearly play an often decisive role in the 
observational studies of co-evolving galaxies and black holes. 
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Fig. 1. — Examples of SED decompositions. Black circles are rest-frame fluxes corresponding to the 
14 bands used to constrain the SED of each object. Purple and blue lines correspond respectively 
to the galaxy and the AGN template found as best fit solution through the minimization for 
the BC03 template set (red and dark green for the P07 one), while the black line shows their sum 
(dark grey for P07 total). Pink and cyan dotted lines show the range of allowed SED from the 
BC03 template library within la of the best-fit Mk measure, and light gray their sum. For one 
objects (VIMOS IDs 834988) our fitting procedure returns only an upper limit for the galaxy rest- 
frame K-band magnitude (overall 10/89); in this case we only plot the AGN spectral component. 
The inset in each panel shows the normalized probability distribution, P = exp (— Xred/^)' 
rest-frame K-band absolute magnitude of the host galaxy, with the solid vertical line marking the 
best fit value and the dashed lines the 1-a uncertainties. The red vertical line marks instead 5% of 
the K-band magnitude of the AGN component. The full set of images of SED decomposition can 



be found at this URL: http://www.mpe.mpg.de/~am/piot_sed_all_rev.pdf 
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Fig. 2. — Estimated rest-frame K band absolute magnitude vs. total stellar mass for the type-1 
AGN hosts of our sample (large panel). Black solid circles are measures, red arrows upper limits. 
As a reference, lines of constant mass to light ratio {M/Lk) equal to 1 and 0.1 are plotted as 
dashed and dot-dashed lines, respectively. The blue line is the average mass to light ratio for 
blue/actively star-forming galaxies at z = 1.5 computed with the Arnouts et al. (2007) relation 
and shifted by -0.1 dex to make it consistent with the S-COSMOS results (Ilbert et al. 2009): 
LogM/Lx = 0.27 z — 0.15. The purple line is the average mass to light ratio for red/quiescent 
galaxies at z = 1.5 computed with the Arnouts et al. (2007) relation, also shifted by -0.1 dex 
(Ilbert et al. 2009): LogM/Lx = 0.17 z — 0.05. The upper small panel shows the distribution of 
Mk (black histogram: detections, red histogram: upper limits) together with the distribution of 
AGN K-band rest frame magnitudes (green histogram). In the right small panel the distribution 
of the measured stellar masses is displayed. 
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Fig. 3. — A few examples of the spectral fitting of the Mgll spectral region. For each object, 
the upper panel shows the fit to the power-law continuum (red dotted line) and the Fe II emission 
(red solid line) that we then subtract to the spectrum. In the bottom panels we show instead the 
same spectrum after the Fe II and continuum subtraction, with the solid green line is the final best 
fit. Note that some of the apparent absorption components redwards of the line are not physical, 
and are required by the fitting routing to compensate for an incorrect representation of the Fell 
emission. In the top left panel, 801709 is one of the three objects in the sample showing signs of 
absorption in the Mgll line region. They are highlighted as blue open circles in Fig. [Sj for them, 
the FWHM is computed from the emission components only. 
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Fig. 4. — Large Top Panel: Black hole mass vs. bolometric luminosity for the zCOSMOS BLAGN 
in the redshift range 1 < z < 2.2. Different colors correspond to different redshift ranges, while 
dashed lines mark the loci of constant Eddington ratios. Large Bottom Panel: Eddington ratio vs. 
bolometric luminosity for the same objects. Dotted lines mark the loci of constant BH mass, while 
the black solid line is the best fit linear regression to the data points, given by Log(Lboi/-^^Edd) = 
—1.38 + 0.64LogLboi,45, where Lboi,45 is the bolometric luminosity in units of 10^^ erg s"^. The 
small panels in the top and right hand side display the distributions of bolometric luminosity, black 
hole mass and Eddington ratio, respectively. 
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Fig. 5. — Scaling relations for zCOSMOS type-1 AGN in the redshift range 1 < z < 2.2. Each row 
show a different redshift interval (lowest 1 < z < 1.25; middle 1.25 < z < 1.6; top 1.6 < z < 2.2). 
Left panels: Black hole mass host K-band absolute magnitude relation. Filled symbols represent 
measurements, leftwards arrows upper limits on the host luminosity. The black solid line is the 
best fit to the Graham (2007) local spheroids sample relation, with dashed lines marking a ±0.3 
dex offset. The typical error bars of our measurements are shown as black cross in the lower right 
corner. Black open circles mark the location of our galaxies when passively evolved down to z = 
assuming a formation redshift of zj = 3. Right panels: Black hole mass host stellar mass relation. 
Open triangles (filled squares) denote the objects with low (< 0.33) and high (> 0.33) Galaxy 
to AGN luminosity ratio in the rest frame K-band, respectively. Symbols with le ftwards arrows 



repres ent upper limits on the host mass. The black solid line is the best fit to the iHaring &: Rix 



(|2004l ) local spheroids sample relation, with dashed lines marking a ±0.3 dex offset. Red squares 
mark the objects detected in the radio by VLA at 1.4 GHz, while red stars mark the location of 
the 5 zCOSMOS AGN in the Jahnke et al. (2009) sample. Blue circles mark the objects with 
absorption features in the Mgll line. The typical error bars of our measurements are shown as 
black cross in the lower right corner. 
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Fig. 6. — Redshift evolution of the offset measured for our type-1 AGN from the local 
Mbh — relation. Different colors identify different ranges of Eddington ratios (purple cir- 
cles Log (-C/boiZ-f'Edd) < — !•) and light blue ones Log (iboi/-^Edd) > — !•) with upwards arrows 
repres enting upper limits on the host mass. The offset is calculated as the distance of each point 
to the iHaring &: Rix I (|2004l ) correlation. Solid black line shows the best fit obtained assuming an 
evolution of the form ALog(MBH/-M*)(-2) = '^2Log(l + z); for which we found 82 = 0.68 ± 0.12. 
The red lines show the bias due to the intrinsic scatter in the scaling relation to be expected even 
if they are universal. Solid line is for an intrinsic scatter of 0.3 dex; dashed of 0.5 dex; dot-dashed 
of 0.7 dex (see text for details). In the inset, we show a comparison of our data (black circles) with 
data from the literature, plotted as green open symbols: triangles are from Salviander et al. (2007, 
low-z) and Shields et al. (2003, high-z); squares from Woo et al. (2008) and circles from Peng et 
al. (2006b). 
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Fig. 7. — Black hole mass host stellar mass relation for zCOSMOS type-1 AGN in the redshift range 
1 < z < 2.2. Symbols with le ftwards arrows r epres ent upper limits on the host mass. The black 
solid line is the best fit to the iHaring fc Rix I (120041 ) local spheroids sample relation, with dashed 
lines marking a ±0.3 dex offset. Red arrows represent the direction of evolution of the points in the 
MBH-Af* plane in 300 Myr on the basis of their instantaneous accretion- and star formation-rates 
and an AGN duty-cycle estimated from the amplitude of the corresponding luminosity and mass 
functions (see text for details). 
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Fig. 8. — Offset from the local A/bh — relation as a function of black hole mass. Filled symbols 
represent measurements, symbols with upwards arrows represent upper limits on the host mass. 
Open triangles (filled squares) denote the objects with low (< 0.33) and high (> 0.33) Galaxy 
to AGN luminosity ratio in the rest frame K-band, respectively. Lines show the bias due to the 
intrinsic scatter to be expected even if the local relation is universal. Solid lines are for an intrinsic 
scatter of 0.3 dex; dashed of 0.5 dex and dot-dashed of 0.7 dex. The typical error bars of our 
measurements are shown as black cross in the upper left corner. 
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Table 1. Rest frame K-band Magnitude and total stellar masses of BLAGN hosts in COSMOS 
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821885 


1 50 708008 


2.292316 


1.099 


-25.78 


0.607 


11.26 







±0U-0c5ioDo 


z.zo / uy ( 




-zo.oy 


u.zoy 


1 1 '^Q 

±1 .oy 


u 


822703 


150.536163 


2.273239 


1.087 


-25.61 


1.854 


10.97 





823199 


150.451859 


2.144812 


1.298 


-25.00 


0.206 


10.80 





823714 


150.345932 


2.147529 


1.258 


-26.16 


0.841 


11.22 





824176 


150.236267 


2.289114 


2.078 


-25.11 


0.147 


10.22 





824306 


150.214676 


2.204261 


1.841 


-24.03 


0.068 


10.08 





824390 


150.199768 


2.190844 


1.510 


-25.79 


0.381 


11.09 





824396 


150.198990 


2.132499 


2.160 


-27.12 


0.223 


11.56 





824572 


150.158371 


2.139555 


1.828 


-26.61 


0.207 


11.29 





825363 


150.004471 


2.237096 


1.407 


-25.77 


0.317 


11.24 





825899 


149.895416 


2.239492 


1.742 


-25.52 


0.179 


10.91 





825906 


149.894852 


2.174454 


1.323 


-25.48 


0.991 


11.13 





827274 


149.624283 


2.180656 


1.185 


-25.93 


0.478 


11.05 





829667 


150.499130 


2.444901 


2.025 


-27.36 


0.560 


11.04 





829682 


150.495667 


2.412547 


1.370 


-26.13 


0.818 


11.18 





830510 


150.347702 


2.390998 


1.848 


-26.02 


0.379 


11.28 
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Table 1 — Continued 



VIMOS ID RA (deg) DEC (deg) z A/^lVega] /gaLK^ Log A/* [A/©] ul*^ 



831077 


150.231812 


2, 


.363971 


1.936 


-25.02 


0.093 


10.53 





832354 


149.993912 


2, 


.301415 


1.789 


-25.56 


0.130 


10.86 





832715 


149.919785 


2, 


.327419 


1.454 


-26.04 


0.136 


11.18 


1 


832923 


149.881012 


2, 


.450839 


1.315 


-26.16 


0.957 


11.14 





832961 


149.871841 


2, 


.342855 


1.735 


-25.19 


0.487 


10.76 





832963 


149.870712 


2, 


.417283 


1.528 


-25.06 


0.353 


10.78 





833273 


149.812943 


2, 


.345459 


1.800 


-26.64 


0.464 


10.67 





833541 


149.763458 


2, 


.334125 


1.131 


-25.24 


0.431 


10.81 





833712 


149.730255 


2, 


.453799 


1.101 


-24.58 


1.592 


10.61 





833817 


149.705872 


2, 


.419752 


1.108 


-26.16 


0.525 


11.17 





834079 


149.660278 


2, 


.410915 


1.161 


-26.28 


0.768 


11.48 





834383 


149.602066 


2, 


.392675 


1.849 


-27.39 


0.215 


11.72 


1 


834988 


149.462845 


2, 


.356840 


1.185 


-24.72 


0.186 


10.65 


1 


835006 


149.459167 


2, 


.430080 


1.242 


-23.80 


0.339 


10.32 





835631 


150.715118 


2, 


.484831 


1.996 


-26.44 


0.117 


10.76 





835840 


150.668900 


2, 


.516766 


1.573 


-23.39 


0.088 


9.72 





836198 


150.597351 


2, 


.617924 


1.447 


-26.12 


1.295 


11.16 





836355 


150.572617 


2, 


.499909 


1.102 


-24.70 


0.190 


10.61 





837652 


150.334427 


2, 


.561485 


1.832 


-26.31 


0.234 


11.19 





837827 


150.305588 


2, 


.602232 


1.342 


-24.98 


0.188 


10.76 


1 


837858 


150.299744 


2, 


.506903 


1.506 


-25.01 


0.222 


10.77 


1 


838223 


150.230820 


2, 


.578165 


1.401 


-25.58 


0.132 


10.99 


1 


838610 


150.163803 


2, 


.597661 


1.589 


-25.09 


0.367 


10.78 





839188 


150.058716 


2, 


.477386 


1.256 


-24.95 


0.095 


10.74 


1 


839751 


149.955612 


2, 


.502021 


1.458 


-25.85 


0.583 


11.24 





841635 


149.575638 


2, 


.575658 


1.171 


-24.50 


0.590 


10.54 





843302 


150.635406 


2, 


.649920 


1.223 


-25.22 


0.343 


10.86 





843685 


150.554871 


2, 


.641009 


1.144 


-25.36 


0.503 


10.81 





844213 


150.456650 


2, 


.648144 


2.050 


-26.43 


0.279 


11.27 





845220 


150.251328 


2, 


.737164 


2.162 


-25.32 


0.104 


10.85 





845272 


150.240799 


2, 


.659021 


1.410 


-24.87 


0.449 


10.86 





845728 


150.147079 


2, 


.717479 


1.177 


-25.62 


1.152 


11.14 





845970 


150.104462 


2, 


.691239 


1.882 


-26.03 


0.272 


11.17 


1 


846335 


150.042480 


2, 


.629174 


1.569 


-26.63 


0.251 


11.46 





847623 


149.774170 


2, 


.674153 


1.108 


-26.02 


1.851 


11.26 





848220 


149.621902 


2, 


.738307 


1.889 


-25.78 


0.088 


10.39 





851007 


150.158981 


2, 


.825123 


1.856 


-26.16 


0.332 


11.23 





900028 


150.084366 


2, 


.290529 


1.112 


-24.79 


0.595 


10.71 





900066 


149.551071 


2, 


.316902 


1.432 


-25.60 


0.462 


11.33 





950021 


149.574463 


2, 


.085072 


1.623 


-25.38 


0.452 


10.86 






^Ratio of galaxy to AGN luminosity in the rest frame K-band, as derived from the SED decom- 
position technique. 

''Upper limit flag: l=upper limit 
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Table 2. Mgll broad emission line FWHM, continuum luminosity and black hole mass for 

zCOSMOS type-1 AGN 



VIMOS ID 




FWHM 1000 


AFWHM 


-^3000, 40*^ 


AZ/3Q00 40 


Log Mbh fMol^ 


Radio Flag'' 


801709 


22.41 


9.23' 


0.19 


2.97 


0.15 


8.67 


RQ 


803695 


21.39 


4.98 


0.06 


8.34 


0.20 


8.35 


RQ 


805949 


21.74 


2.52 


0.15 


5.39 


0.21 


7.67 


RL 


807560 


20.36 


7.82 


0.21 


37.48 


0.48 


9.05 


RQ 


808150 


20.20 


3.45 


0.11 


85.33 


5.14 


8.51 


RQ 


810061 


21.71 


4.69 


0.41 


18.64 


1.10 


8.46 


RQ 


811239 


20.83 


2.63 


0.05 


33.89 


0.92 


8.08 


RQ 


811646 


21.99 


5.48 


0.41 


15.78 


1.51 


8.56 


RL 


811799 


22.41 


3.80 


0.29 


9.21 


0.37 


8.13 


RQ 


811904 


22.22 


3.74 


0.56 


2.89 


0.18 


7.88 


RQ 


811960 


20.90 


5.46 


0.81 


33.60 


0.97 


8.71 


RQ 


813283 


20.88 


4.97 


0.21 


21.60 


0.30 


8.54 


RQ 


813416 


22.42 


3.63 


0.25 


5.98 


0.56 


8.01 


RQ 


813886 


22.42 


6.49 


0.28 


3.18 


0.21 


8.38 


RQ 


814414 


20.34 


2.70 


0.12 


26.15 


0.39 


8.05 


RQ 


816818 


21.14 


5.46 


0.07 


11.64 


0.27 


8.50 


RQ 


817202 


21.27 


3.80 


0.10 


18.89 


0.54 


8.28 


RQ 


817260 


20.95 


6.77 


0.09 


45.04 


1.80 


8.96 


RQ 


817480 


22.30 


10.46 


0.26 


11.40 


1.59 


9.06 


RQ 


818094 


20.13 


4.42 


0.17 


87.15 


2.41 


8.72 


RL 


819187 


21.03 


5.29 


0.06 


33.82 


1.63 


8.69 


RQ 


819193 


19.00 


4.46 


0.21 


221.60 


2.12 


8.92 


RQ 


819446 


21.73 


3.35 


0.15 


14.75 


1.53 


8.12 


RQ 


819579 


21.69 


3.95 


0.14 


7.83 


0.27 


8.13 


RQ 


819592 


21.90 


3.82 


0.11 


6.55 


0.21 


8.07 


RQ 


819644 


18.53 


3.65 


0.45 


139.86 


6.38 


8.65 


RQ 


819702 


20.93 


3.56 


0.16 


24.83 


0.63 


8.28 


RQ 


820341 


20.62 


4.69 


0.58 


22.64 


0.31 


8.50 


RQ 


820375 


22.45 


5.02 


0.16 


8.07 


0.77 


8.35 


RQ 


820673 


21.06 


4.97 


0.07 


25.19 


1.34 


8.57 


RQ 


820679 


19.82 


5.39 


0.07 


63.30 


0.65 


8.83 


RQ 


821039 


22.01 


5.97 


0.43 


4.79 


0.20 


8.39 


RQ 


821885 


20.99 


5.48 


0.19 


12.85 


0.22 


8.52 


RQ 


822461 


21.29 


7.93 


0.52 


14.01 


0.32 


8.86 


RQ 


822703 


21.58 


5.88 


0.31 


6.60 


0.22 


8.45 


RQ 


823199 


20.61 


3.94 


0.13 


26.23 


0.40 


8.38 


RQ 


823714 


20.95 


2.68 


0.21 


19.24 


0.34 


7.98 


RL 


824176 


21.48 


6.19 


0.19 


25.54 


2.31 


8.77 


RL 


824306 


20.99 


2.68 


0.10 


27.69 


1.14 


8.05 


RQ 


824390 


21.12 


4.90 


0.27 


19.09 


0.65 


8.50 


RQ 


824396 


19.49 


3.62 


0.19 


108.25 


4.59 


8.59 


RQ 


824572 


20.45 


4.27 


0.47 


52.42 


1.44 


8.59 


RL 


825363 


21.70 


4.54 


0.23 


10.70 


0.36 


8.32 


RQ 


825899 


20.85 


3.50 


0.92 


32.38 


1.96 


8.32 


RQ 


825906 


22.05 


6.09 


0.10 


6.52 


0.24 


8.47 


RQ 


827274 


20.82 


2.89 


0.05 


17.00 


0.37 


8.02 


RQ 


829667 


19.29 


5.60 


0.27 


218.38 


3.40 


9.12 


RQ 


829682 


21.60 


4.97 


0.04 


12.63 


0.21 


8.43 


RQ 


830510 


22.17 


6.76 


0.23 


17.24 


0.85 


8.76 


RQ 
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Table 2 — Continued 



VIMOS ID 




FWHMiooo'' 


AFWHM 


3000, 40*^ 


AL3ooo,40 


Log Mbh [MQ]d 


Radio Flag<= 


831077 


20.58 


4.13 


0.11 


55.50 


3.90 


8.57 


RQ 


832354 


20.12 


4.45 


0.09 


67.77 


1.49 


8.68 


RQ 


832715 


20.15 


7.90 


0.14 


50.76 


0.81 


9.12 


RL 


832923 


21.29 


3.59 


0.26 


15.58 


0.35 


8.19 


RL 


832961 


21.86 


3.33 


0.24 


12.29 


0.71 


8.08 


RQ 


832963 


22.40 


4.93 


1.25 


6.90 


0.55 


8.30 


RQ 


833273 


21.94 


5.23 


0.38 


14.22 


1.83 


8.50 


RQ 


833541 


21.11 


4.99 


0.30 


10.49 


0.34 


8.40 


RQ 


833712 


22.03 


4.76 


0.22 


2.43 


0.44 


8.06 


RQ 


833817 


20.57 


5.56 


0.39 


27.18 


0.45 


8.69 


RQ 


834079 


21.41 


11.18 


0.53 


7.70 


0.31 


9.03 


RL 


834383 


19.86 


5.32 


0.12 


111.41 


1.68 


8.94 


RQ 


834988 


20.57 


5.53 


0.21 


18.92 


0.39 


8.61 


RQ 


835006 


22.20 


1.93 


0.10 


5.64 


0.28 


7.45 


RQ 


835631 


19.54 


2.31 


0.10 


144.30 


2.77 


8.26 


RQ 


835840 


22.29 


3.19 


0.33 


6.17 


0.75 


7.90 


RQ 


836198 


21.76 


3.92 


0.12 


10.76 


0.35 


8.19 


RQ 


836355 


20.68 


5.10 


0.14 


17.32 


0.22 


8.52 


RQ 


837652 


20.08 


5.90 


0.11 


76.87 


1.86 


8.95 


RQ 


837827 


22.06 


3.82 


0.10 


6.02 


0.28 


8.05 


RQ 


837858 


21.02 


3.54 


0.09 


24.02 


1.40 


8.27 


RQ 


838223 


20.05 


4.31 


0.06 


54.27 


0.92 


8.61 


RL 


838610 


21.88 


3.94 


1.10 


12.12 


1.13 


8.22 


RQ 


839188 


20.55 


4.83 


0.07 


28.03 


0.36 


8.57 


RL 


839751 


22.11 


9.04 


0.74 


5.90 


0.76 


8.80 


RQ 


841635 


21.72 


5.80^ 


0.41 


7.01 


0.23 


8.45 


RQ 




21 15 


o.oz 


21 


1 fi 79 


n 97 


R 9fi 




843685 


21.53 


5.55 


0.11 


6.48 


0.22 


8.39 


RQ 


844213 


20.47 


4.76 


0.36 


62.15 


2.02 


8.72 


RQ 


845220 


20.85 


3.92 


0.20 


25.43 


3.22 


8.37 


RQ 


845272 


22.49 


3.79 


0.26 


5.59 


0.42 


8.03 


RQ 


845728 


22.19 


6.19 


0.16 


3.43 


0.29 


8.36 


RQ 


845970 


21.08 


3.79 


0.20 


35.96 


2.41 


8.41 


RQ 


846335 


20.21 


4.30^ 


0.28 


57.32 


1.29 


8.61 


RQ 


847623 


21.59 


6.54 


1.44 


4.55 


0.23 


8.46 


RQ 


848220 


20.03 


3.83 


0.19 


94.11 


2.76 


8.62 


RQ 


851007 


20.72 


3.38 


0.11 


41.31 


1.23 


8.34 


RQ 


900028 


22.49 


8.21 


0.24 


1.66 


0.20 


8.45 


RQ 


900066 


22.36 


6.89 


0.60 


6.14 


0.25 


8.57 


RQ 


950021 


22.11 


6.43 


0.17 


12.09 


1.77 


8.65 


RQ 



''Apparent i-band (AB) magnitude of the zCOSMOS spectroscopy target 
•'Mgll emission line FWHM in units of 1000 km s~l 
'^Monochromatic continuum AGN luminosity in units of 10*" erg s~^ 

'^Logarithm of the black hole mass, as computed with Eq.([TJ, adopting the MG08 calibration (see text for 
details) 

<=RQ are objects detected in the VLA/COSMOS Survey at 1.4 GHz llBondi et al. IbOOSi : [Schinnerer et al. Il2007h 
'Objects with clear narrow absorption features in the Mgll line 



